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1 Introduction 



Compact Einstein spaces with positive-definite metrics have been studied extensively: 
first, the spaces are considered as candidates of compact internal spaces which are admit- 
ted in higher-dimensional gravitational theories [I]; and secondly, the spaces are expected 
that they may dominate the path integrals of quantum gravity 0. 

There are many examples of homogeneous Einstein metrics to be found in a lot of 
literature, but they are very exceptional among general Einstein spaces. In contrast, the 
inhomogeneous Einstein metrics would cover wide range, but our knowledge for them 
seems to be largely limited. The first example of inhomogeneous compact Einstein space 
which is a solution to the Einstein equation with a positive cosmological constant, 



was constructed by Page[3j by taking a limit of the Euclidean Kerr-de Sitter solution, 
and it was generalized by Berard-Bergery[3]. 

In the absence of any general understanding of the solution to the Einstein equa- 
tion (Jl.lj) . one of the standard strategies for constructing inhomogeneous examples is to 
study the space with cohomogeneity one metric: the space admits a Lie group action by 
isometries whose orbits span the space with codimension one. It is considered that the 
space which foliates into a sequence of homogeneous subspaces of codimension one. The 
Einstein equation for such metric reduces to a set of non-linear second order ordinary 
differential equations[5j. 

It is possible to replace the homogeneous subspaces by spaces with bundle structure. 
This idea originates from the Kaluza-Klein construction. The Einstein equation trans- 
lates into a coupled system of equations involving the Ricci curvatures of the fibers and 
base space, as well as the curvature of the connection. When we have a suitable choice 
for the bundle space and connection, the Einstein equation also reduces to a set of the 
ordinary differential equations. 

In this paper, we construct compact inhomogeneous examples of the Einstein met- 
ric with a positive cosmological constant in 4n+4 dimensions on the spaces with bun- 



Rab — A<7o6 
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die structure. More precisely we consider the union of principal S'0(3)-bundles P over 
Quaternionic Kahler manifolds, so our spaces may be identified locally with I x P for 
some interval I of the real line. This geometrical setting was studied in jH] and they 
constructed a family of inhomogeneous compact Einstein spaces. Our construction is a 
generalization of their system |S]. 

In this framework, it is possible to take three types of boundary conditions at two 
endpoints of / for completeness, i.e., two types of bolt singularities associated with the 
Quaternionic Kahler manifold and with its twistor space respectively, and the nut singu- 
larity describing S' 4n+3 -collapsing. We find new compact Einstein spaces which have the 
two types of bolts, numerically. The new solutions together with the known solutions give 
a unified description of compact Einstein spaces with a positive cosmo logical constant. 

The paper is organized as follows: Section 2 contains our Kaluza-Klein metric ansatz 
and calculations of the Riemannian curvature. In section 3 we derive the Einstein equa- 
tion after a short review of Quaternionic Kahler manifolds. We also prove the diago- 
nalizability of the metric using the technique for the 4-dimensional Bianchi IX metric. 
In section 4, we discuss the boundary conditions and give asymptotic solutions near 
the boundaries. In section 5, known exact solutions are listed. In section 6 we present 
new solutions obtained by numerical integrations. Section 7 is devoted to summary and 
discussion. 

2 Metric Ansatz and Curvature Calculation 

In this section, we shall consider metrics on (m + 4)-dimensional manifolds M. To make 
the analysis manageable, we assume the following geometrical condition for M. Let 
7i : P — > M be a principal >SO(3)-bundle over an m-dimensional Riemannian manifold 
(M,qm) and be an S'0(3)-connection on P. The connection locally takes the form 

= s^As + s^ds, seSO(3). (2.1) 
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Here, A is an so(3)-valued local 1-form on M and s _1 c?s is considered as the Maurer- 
Cartan form. Let <p l be the components of for the standard basis {Ei} of so(3) which 
satisfies the Lie bracket relations [Ei,Ej] = Yl\=i e ijkEk- By using the left-invariant 
1-forms <7 l defined by s~~ l ds = Y^=i ^Ei, the equation (|2.1j) can be also written as 

3 

^ = Y,A j jt + a l (2.2) 

3=1 

with the adjoint representation 

3 

S^EiS = & i3 E 3 ■ ( 2 - 3 ) 
3=1 

In this setting we consider metrics on M which is locally the product space I x P, where 
/ denotes some interval of R. Given a metric 6^- on SO(3), the Kaluza-Klein metric takes 
the form 

3 

g = dt 2 + + f(t) 2 9M ■ (2.4) 

i,3=l 

We can show that the matrix bij(t) is diagonalizable for all t under a certain condition 
of the base space M (see a proposition in section 3). Thus we write the metric as 

g = dt 2 + o(f) W + b(tf{<p 2 ) 2 + c(t) 2 (0 3 ) 2 + f{tfg M . (2.5) 

If we impose the condition a = b = c, then the metric ()2.5|) has an isometry 5*0(3). 
Otherwise it has no such symmetry because of the explicit dependence on the group 
element 0y. 

Now, we choose an orthonormal basis e a = {e°, e l , e a ; % = 1 ~ 3, a = 1 ~ m} for 
(ESI), 

e° = dt, e l = (nfi, e a = fe a , (2.6) 
where (cij) = (a, 6, c) and e a is an orthonormal basis for g^- Then, the spin connection 
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uj a b defined by de a = 



—ub a b A e b is calculated as 



./ 



1 

a i J & 



a\ - aj - a 2 fc 

/ j £ ijk _ _ e ? 



(2.7) 



a,- 



J k,/3 k,i •* 



where F = J2 i F' l Ei is the curvature 2-form of A and Lu a/3 is the spin connection on M. 
The curvature 2-form fl a b = di2j a b + J2 C ^ ac ^ <^ cb * s §i ven by the following equations: 



^ «i „ J- 



(Xj GLj Otj^ CL^ CLj 



— + 

CLjQ-k CLj CLiCLjCLfc 



& A e k 



ai cii f 



V 2 a t f 



k,a,0 



f 

a =--e°Ae°- 



&i I d{ f 



f 



2/ 2 a t f 



k,i,f3 

^7 ^' 5a/3 + Z^4p £ ^ r 



a,/ ^ 1 af + a 2 - a 2 



a-if 

i 

2 . /2-™ 7 



k,l,-y 



e j A e' 3 



k,p 1 Kp n 
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f 2 a 2 
^ = ^ -yA^-J] {®^F k aP @uF\ s + ki F^0 H F l ps ) e< A e 5 

' i,k,l 7,<5 

i,k 7 »,fc V / 



jj fc,Z 7 



1 



2/ : 



5>i 



fl2 + fl 2 _ fl 2 2 + fl 2 _ a 2 

0ki e A e + fc2 



■e'^Ae 1 



a 2 a 3 



a 3 ai 



+0fc 3 e 1 A e 2 

a^ 



where Qq,^ is the curvature on M and D 1 FK is the gauge covariant derivative of F, 



1,3 



The remaining components Qij = 1 ~ 3) are 



(2.8) 



n 
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a 3 ai a\ + a 2 — a 2 , a 2 a\ — a\ + a 2 . 



a\a 2 ai 2aia 2 a3 a 2 2a\a 2 a 3 



dia 2 a 2 — a\ — a 2 a\ — (a\ — a 2 ,) 2 
+ + 



a,\a 2 



2ala 2 , 



e°Ae 3 



e'Ae 2 



as — a? — as 



n ~ "2 fc a l fl 2 >r^„ fe „ / 

— — Wk3* a p SP^n^nfi ^ 

a l a 2 J , 

',7 



(2.9) 



e a A e/3 , 



and Q 2 3, ^31 ar e given by cyclic permutations 1 — > 2 — > 3 — > 1. The non-zero components 
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of Ricci tensor R a b become 

— - . — , ri - 

Roo 



m 

n- 



f 



R 





di 1 


\ai) 





dk 



^•+4 O(3) + ^EE ^A' 



k,l a,7 



R, 



a/3 



- ~ (m - l) f j 2 - J2 ^ ) Saf, + J 2 Rap ~ E E £i ^U F k > 

' " i,k,l 7 



f 



P 



2f 



R 



a,i 



2/ 3 



k 

ap- 



k,P 



(2.10) 

SO(3) 

Here i? Q/ g denotes the Ricci tensor on M and R^ the Ricci tensor on 5*0(3). Explicitly 
the non-zero components of R^°^ are given by 



a\ — (a 2 , — a^ 2 



2„2„2 ' 22 



at — (a 2 . — a 2 ^ 2 



2„2„2 ' ""-33 



R% 



at — (a 2 — a^ 2 



SO (3) _ 1 "3J r>SO(3) _ 2 V"3 "17 p 50(3) _ 3 V"l "2 



2afa2ag 



, (2.11) 



3 Einstein equation based on Quaternionic Kahler 
Manifold 

In order to solve the Einstein equation, we need a further assumption on the base space 
M of the principal bundle P. Then we will obtain a generalization of the ordinary 
differential equations studied in [HI IE] • 

Let (M,qm) be a 4n-dimensional Quaternionic Kahler manifold (QK manifold). It 
has a set of three almost complex structures J a (a = 1 ~ 3) which satisfy the quaternion 
algebra 

J a J b = -5 ab + e abc J c , (3.1) 
and one can find local 1-forms A 1 such that 



v r = e abc A b ® r 



(3.2) 
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where V denotes the Levi-Civita connection [HJ. According to ^U], we introduce an or- 
thonormal basis with 2-indices {e Ml }; 

3 n 

Then the 2-forms J a are given by 

3 

i=l 
3 

J 2 = ^ e 0i Ae 2l + e 3i A e H , (3.4) 

8=1 

3 

J 3 = e 0i A e 3i + e u A e 2i . 

8=1 

Furthermore, the 5'0(3)-connection A = $^f =1 A*.Ei is Ci-self-dual in the sense of [TU] . 
i.e., the Yang-Mills curvature F = dA + A A A satisfies the following self-dual equation 

3 

*F = ciFAO"- 1 , Q = J2 jaAja ( 3 - 5 ) 

a=l 

with Ci = 6n/(2n + 1)! . We note that the connection A automatically satisfies the 
Yang-Mills equation like a 4-dimenional ordinary instanton. In fact 

D * F = Cl D(F A W 1 - 1 ) = , (3.6) 

where we have used dVt = by ()3.2j) and the Bianchi identity DF = 0. It is known that 
any QK manifold is Einstein and F simply takes the formjnjjH] 

F' = ^- 2 J\ (3.7) 

where A is the Einstein constant for the metric gu- 

Now let us turn to the evaluation of the Ricci tensor (j2.1()jl . From QK geometry, we 
saw that the S'0(3)-connection A is Ci-self-dual and its curvature satisfies the quaternionic 
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relations. Thus, when we take a 4n- dimensional QK manifold as the base space with self- 
dual connection, we have 1 

R al3 = (n + 2) 6 Qf} , J2 D * F ^ = ( 3 - 8 ) 

a 

and 



a ' P (3.9) 

k,l 7 k,l 

These equations finally lead us to the (4n+4)-dimensional Einstein equation with a cos- 
mological constant A: 

a b c f 

-+ 7+ - + 4n-= -A . 
a o c j 

a a lb c A a 4 - (b 2 - c 2 ) 2 a 2 

a = -a^ + c + 4n /J + 2a W + "7^ A ' 

b b (c a A 6 4 - (c 2 - a 2 ) 2 6 2 

6 = " 6 + a + 4n j) + 2aW + " A ' (3-10) 

c c /a 6 A c 4 - (a 2 - 6 2 ) 2 c 2 

c = " a + ^ + 4n 7 + 2a W + "T 4 " A ' 



/ //a 6 c f\ a 2 + b 2 + c 2 n + 2 

7=-7U + 5 + ^ +(te - 1) 7 ~ + ~ 



Here, if we impose the conditions a = b = c ^ f and a ^ b = c = f , then (j3.1()|) yields the 
equation studied in [S] and [7j, respectively (see section 5). Also the singular reduction 
a = and b = c gives the (4n + 3)-dimensional Einstein equation discussed in jH]. In 
case of the trivial bundle P = M x 50(3), i.e., A = 0, the Ricci tensor (j2.10)) directly 
yields the Einstein equation without the assumption of QK manifolds. This equation is 



1 We have used a normalization A = n + 2 in (SIB) an d (|3.9|) since we shall consider only A > cases. 



s 



simply given by dropping the / 4 terms and replacing the dimension An by arbitrary one 
in (EHUj) . 

It is worth noting that the 8-dimensional Einstein equation with vanishing A is spe- 
cial in the sense that the Spin(7) holonomy condition leads to the following first-order 
equations [I2J Uni; 



a 


a 2 -(b- 


cf 


a 


a 


2abc 






b 


b 2 -(c- 


a) 2 


b 


b = 


2abc 
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c 


c 2 — (a — 


b) 2 


c 


c 


2abc 






f 


a + b + c 






f~ 


2f 2 ■ 







(3.11) 



We can verify by substituting ()3.11|) into ()3.10|) with n = 1 that the metric is indeed 
Ricci-flat, and the explicit solutions were constructed in [SJ [T3J EH EC3 • 
Finally we prove the following as mentioned in section 2 : 

Proposition Let (M, g^) be a 4n-dimensional QK manifold with ci-self-dual 
connection. Then the Kaluza-Klein metric (|2.4|) for the Einstein space can be put in the 
diagonal form (J2.5|) for all t. 

Proof. We calculate the Ricci tensor for the metric ()2.4|) . and the result is given in 
Appendix A. Note that the off-diagonal component i?o; takes the same form as the 
4-dimensional Bianchi IX type cosmological model. So this enables us in the usual way 
(see PHI, for example) to diagonalize the metric for all t . The argument is based on an 
important property of the Einstein equation, namely invariance under the right action 
of 5*0 (3). Indeed, using the corresponding transformation of the connection 

0-^0 = so V so , s E S0(3) , (3.12) 
we can diagonalize the fiber metric 6^ at an initial time t = to. Then the Einstein 
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equation R 0i = leads to fry = for i ^ j at t = to, which implies the diagonality of the 
solution for all time. 



4 Boundary Condition 

In this section, we discuss boundary conditions for the Einstein equation Let us 

assume the following compact conditions of M ~ I x P : 

1. I is the closed interval [ti,t2] , 

2. QK manifolds have a positive scalar curvature, namely A > 0. 

Furthermore we require that the singularities at the boundaries t\ and t% are resolved by 
bolts or nuts; there are three types of resolutions, nut «-> nut, bolt <-> nut and bolt <-> 
bolt (see Fig.l). This means that near the boundary P is locally of the form 



where the radius of round /c-sphere S k tends to zero at the boundary and the ^-dimensional 
manifold B e remains non-vanishing. We find there are three choices of the manifold B l 
consistent with the Einstein equation: (Bl) QK manifold M (£=4n), (B2) twistor space 
Z of the QK manifold (£=4n+2), (B3) empty. In the case of (Bl) or (B2) the singularity 
can be resolved by bolt, we call these singularities QK-bolt and T-bolt, and the case (B3) 
by nut. We summarize these boundary conditions as follows : 



P — > S k x B e (k + £ = 4n + 3) , 



(4.1) 



(Bl) QK-bolt 



$ _> dt 2 + t 2 ((0 x ) 2 + (0 2 ) 2 + (0 3 ) 2 )/4 + a 2 g M , f - . 



(4.2) 



Qm denotes a metric on a QK manifold M. 



(B2) T-bolt 



g^dt 2 + kh 2 ^ 1 ) 2 + p 2 g z , t -> . 



(4.3) 
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gz denotes a metric on a twistor space Z , 

g z = (0 2 ) 2 + (0 3 ) 2 + «V , (4.4) 

which shows Z is an S^-bundle over M, and the metric is Kahler-Einstein for a — 1. 



(B3) nut 



£ _> ^ + ^ ((0 1)2 + (02) 2 + (03) 2 + ^ )/4; t ^ Q _ 



(4.5) 



Here we have translated the boundary t — t\ or t 2 to the origin t — 0, and ct , /3 and are 
free parameters. In the case (B3) the QK manifold is required to be HP(n), otherwise 
it yields the curvature singularity at the boundary t — since it does not describe the 
S' 4n+3 -collapsing (nut singularity). On the other hand, in the cases (Bl) and (B2) the 
manifold M would have an orbifold singularity at t — though arbitrary QK manifolds 
are allowed. Indeed, ^j(0*) 2 /4 represents the metric on 5*0(3) ~ S* 3 / 'Zi rather than 
SU(2) ~ S* 3 , and the range of the angle Q = kip ((f) 1 = dip for the fixed twistor space 
coordinate) does not mean < < 2n generally. For the case (Bl), if we choose the QK 
manifold HP(n) , then the orbifold singularity disappears by lifting to an S , f/(2)-bundle 
(see section 5) . 

Using the Einstein equation we find the following asymptotic behavior of the metric 
near the boundary: 



(Bl) QK-bolt 



a = t/2 + a 3 /6 t 3 + • • • , 
b = t/2 + b 3 /Q t 3 + • • • , 
c = t/2 + c 3 /6t 3 + ■■■ , 
/ = fo + / 2 /2 t 2 + ■ ■ ■ . 



(4.6) 



Here three of /o , 03 , 63 and C3 are free parameters which satisfy the relation 



n(2 + n) 1 
a 3 + fr 3 + C3+ 2 +-(l-w)A = , 



(4.7) 







and the remaining coefficients are determined by them. 
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p 




Figure 1: Compact manifolds M ~ I x P. I is [ti,^] and P is a principal 50(3)- 
bundle over a Quaternionic Kahler manifold. Endpoints at t\ and t 2 are boundaries 
which can take nut or bolts. 



(B2) T-bolt 



a = kt + a 3 /6 t 3 H 

b = c = b + b 2 /2 t 2 + 
f = fo + f 2 /2 t 2 + • • • 



(4. 



where k , bo , fo are free parameters and the remaining coefficients are determined 
by them. The cases k = 2 and 1/N for a positive integer N are exceptional and 
the expansion above must be modified by the following equations : 

1. k = 2 

a = 2t + a 3 /6 t 3 H , 

6 = fe + 6i * + & 2 /2 t 2 + • • • , 

(4.9) 

c = b Q -b l t + b 2 /2 t 2 + • ■ ■ , 

/ = /o + / 2 /2t 2 + --- , 
where &o , &i > /o are free parameters and remaining coefficients are determined by 
these. 
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2. k = 1/N 



a = t/N + a 3 /6 t 3 + ■ ■ ■ , 
b = b + 6 2 /2 t 2 + ■ ■ ■ , 




2iV 



+ ••• 



(4.10) 



/ = /o + / 2 /2t 2 + -- , 



where b , f , 8^ = b 2 N — c 2 n are free parameters and remaining coefficients are 
determined by these. Derivation of (|4.1()jl is viewed in Appendix B. 



c = t/2 + C3/6 t 3 + ■ • • , 
/ = t/2 + f 3 /6 t 3 + ■ ■ ■ , 
where a 3 , 63, c 3 and / 3 are parameters satisfying a 3 + b 3 + c 3 + 4n/ 3 + A/2 = . 

5 Example 

The quaternionic projective space HP(n) is a typical example of QK manifolds. The 
Hopf fibration 5* 4n+3 — > HP(n) is a principal S*?7(2)-bundle over HP(n). It has a 
natural connection such that its horizontal space is the orthogonal complement to the 
fiber with respect to the standard metric on S 4n+3 . This connection is Ci-self-dual[17 
and in case of n = 1 the connection is the well-known BPS instanton. 

We explicitly calculate the Kaluza-Klein metric based on HP(ra) according to sections 
121 and E3 For the base space M=HP(n), the standard metric is written as 



(B3) nut 



a 



b 



t/2 + a 3 /6t 3 + ■■■ , 
t/2 + b 3 /6 t 3 + ■ ■ ■ , 



(4.11) 



9m = 



Adx A dx A Ax A dx A x B dx J 
1 + x c x c (1 + x c x c ) 2 



(5.1) 
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where x A = x A + xfi + x A j + x A k (A = 1 ~ n) are quaternionic coordinates and x A are 
their conjugates. Then the Ci-self-dual connection takes the form 

AH + Aij + A*k= X dX l + J X xB X ■ (5.2) 

Thus we have a Kaluza-Klein metric (j2.5|) on / x P explicitly, and the Einstein equation 
is given by ()3.10|) . 

By the construction P is an S'0(3)-bundle over HF(n). In general, there is an 
obstruction to lifting an 5'0(3)-bundle over QK manifold to an S*?7(2)-bundle, i.e., the 
Marchiafava-Romani class e. A result of jTl] says that e = if and only if the QK manifold 
is HP(n), and so in this case the total space P lifts to the covering space P as an SU(2)- 
bundle. Actually P is the total space of the Hopf fibration, S 4n+3 ~ Spin + 1)/ Spin). 

For completeness we present a summary of the known compact complete metrics based 
on the Hopf fibration. All these metrics are given as solutions to appropriate reductions 
of our equation. 

(Rl) reduction to one variable ; a — b — c — f 

a(t) = 1/2 cost , -7r/2<t<7r/2 (5.3) 

with A = An + 3 . This solution, which obeys the boundary condition (J4.5|) at 
t = ±7r/2 , represents the standard metric on M = S 4n+4 . 

(R2) reduction to two variables ; a = b = c and / 
An explicit solution is given by 

a(t) = 1/2 sin t cost , f(t) = 1/2 cost , < t < vr/2 (5.4) 

with A = 4n+12 . This solution, which obeys ()4.2j) . ()4.5|) at t = 0, 7r/2, respectively, 
gives the standard metric on HP(n + 1) . Another solution satisfying ()4.2j) at the 
both endpoints was constructed numerically. It gives a metric on an S" 4 -bundle over 
HP(n), i.e., M = UP(n+ 1) (j HP(n+ 1)0. The case n = 1 has been analytically 
established by B6hm[~ 



14 



(R2) reduction to two variables ; a and b = c = f 

The situation is very parallel to (R2), but the topology is different as we will see 
shortly. The Fubini-Study metric on CP(2n + 2) belongs to this class. Explicitly, 
the metric is written as 

a(t) = 1/2 sin t cost , b(t) = 1/2 cost , < t < vr/2 (5.5) 

with A = 4n + 6 . This solution obeys the boundary conditions ([4.3)1 . (|4.5jl at 
t = ,7r/2 , respectively. In this case the general solution was constructed by 
Page- Pope |7j as 

ft " (1 " r2)2 " +1 rfr 2 + t^^W) 2 + m^l - r 2 )^ z , (5.6) 



P(r) (i _ r i\ 

where 



P = m2r + ^±lg 2n+1 - A Q 2n+2 , Q n = V f n ) (5.7) 



with integration constants (i=l ,2) . Its special case satisfying (|4.3|) at the both 
endpoints gives a metric on an S 2 -bundle over the twistor space Z, i.e., M = 



CP(2n + 2) t) CP(2n + 2) 



6 New Solutions 

The solutions listed in the previous section have two of three possible boundaries at 
the endpoints: (Bl) QK-bolt, (B2) T-bolt and (B3) nut. Correspondence between the 
solutions and boundaries is schematically shown in Fig.2. No solution is found which 
connects QK-bolt and T-bolt though the existence is naturally expected from the Fig. 2. 
The solutions with QK-bolt are found in the ansatz (R2), and the solutions with T-bolt 
are in the other ansatz (R2), then we search new solutions which connect QK-bolt and 
T-bolt under a generalized ansatz in which the metric has three unknown variables. 
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nut 




QK-bolt 



T-bolt 



Figure 2: Relation between boundaries and solutions. 

We assume b = c, hereafter, and treat a, b, f as unknown metric functions of t. From 
(|3.1U|) we obtain the evolution equations 



u alb f \ a a , 



bab f\ a 2 -2b 2 b 2 
-b[a + b +An j)--2b^ +n T*- A > 



f a b , A a 2 + 2b 2 n + 2 

f 1 \ a + b + ~ / J 2f~ A ' ~p 

Taking a combination of (j3.10|) we also get the Hamiltonian constraint 



(6.1) 



b I a b \ f I a b , f 

2 l [ 2 a + l) +inL ; [ 2 -a +i l + {in - 1)L f 

a 2 -4b 2 a 2 + 2b 2 A . n . 1 ' , A ^ , 
— — + n — — An{n + 2) — + (4n + 2) A = 0. 



2b A 



(6.2) 



We should integrate the equations ()6.1|) from QK-bolt at t = t\ = as an 'initial' 
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value problem in t. Since a and b are zero at t = 0, the set of ordinary differential 
equations is singular there. To deal with the singularity, we use the asymptotic 

form of the solutions 1)4.6 j] near the singularity, and start numerical integrations at an 
initial time t = ti = + e, where e is a small amount of time duration. 

Equation ()4.6|) gives the asymptotic behavior of a, b and / near QK-bolt: 

a = t/2 + a 3 /6 t 3 H , 

b = t/2 + 6 3 /6 t 3 + ■ ■ ■ , (6.3) 
/ = / + / 2 /2t 2 + - , 



where f , a 3 and 63 satisfy 



n(2 + n) 1 



a 3 + 2& 3 + —r 2 +-(l-w)A = 0, (6.4) 



and /2 is given by 



4 V /« 

Introducing a parameter ^ by 



h = U^-Vo)- (6-5) 



^o = ^(&3-a 3 ), (6.6) 



we can solve (|6.4j) as 



1 / 1 n , . 1 
a 3 = -{n - 1)A - -{n + 2 W - 2^ , 
6 6 / 2 

1 / w n , \ 1 
63 = - n - 1 A - -{n + + ^o- 
6 6 /(f 



(6.7) 



Therefore, the initial condition at t — t j is given by (|6.3|) (|6.5|) and (|6.7j) with t = e. We 
have two free parameters /o and ^0 to specify the initial condition near QK-bolt. 

We are interested in finding solutions for which b and / stay positive and finite, and a 
returns to zero at some moment t = ti- At this point the set of equations (|6.1|) becomes 
singular again. We assume that this singularity is resolved by T-bolt ()4.3j) . If k ^ 1 the 
cone-type singularity appears at t — ti so we set k = 1 for regularity. Then the condition 
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()4.10|) with b = c (5\ = 0) shows that the asymptotic behavior of a, b and / near t = t 2 
should be 



a=(t 2 -t) +a 3 /6 (t 2 -t) 3 + --- , 
b = b + b 2 /2 (t 2 -t) 2 + --- , 
/ = /o + / 2 /2 {t 2 -tf + --- , 



(6.8) 
(6.9) 
(6.10) 



where bo and fo are free constants and 03, b 2 and f 2 are described by them as 




(6.11) 



We stop numerical integrations at tf = t 2 — e (e is a small constant) and check whether 



A solution for (|6.1|) defined on the region [tj, tf] gives a map {(fo,9o)} {(V 1 ^ 2 )}- 
Then we can regard V = (V 1 , V 2 ) as a vector field on a 2-dimensional plane parameterized 
by (fo,go)- If V — 0, a and b have the forms of (|6.8|) and (|6.9|) near t = t 2 , in addition, 
it is shown from (|6.2j) / should have the form of (J6.1(J|) automatically. Then, vanishing 
of V means the endpoint at t = t 2 is T-bolt. 

Numerical integrations are done by using a fourth-order Runge-Kutta routine. We 
have verified that the constraint ([6.2)1 . which is preserved by the evolution equations 
(|6.1|) . holds in high accuracy in our numerical integrations. We also reproduce known 
solutions listed in the previous section when we set a = b or b = f. 



the values of a, b, f and a, b, f agree with ()6.8|) - ()6.1()|) . 
Here, we define 



V 1 = a(t f ) - a T , 
V 2 = b(t f )~b T , 



(6.12) 



where 



clt = — 1 — ^3/2 e 2 , bx 



b 2 e. 



(6.13) 
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5 10 15 20 

trtf 

Figure 3: Vector field V on (/ , go)-plane forn = 1 case. Parameters mf and m g give 
fo and g as f = 1. 75333405+1. Ox xm/ and g = 0.002532820+1. Ox 10~ 9 xm g . 
Horizontal and vertical components of the arrows are V 1 and V 2 , respectively. The 
vector field vanishes on a point (/ *,^) ~ (1.75334014,0.002532828). 

Fig. 3 shows the vector field V on (fo, go) -plane for n = 1 case. We normalize A = 1. 
There is a critical point (fo,9o) on which V vanishes. The map (/o,<7o) ^ V = 
corresponds to the solution for ()6.1|) which connects QK-bolt at t = and T-bolt at 
t = t 2 - Evolution of a, b and / in t is plotted in Fig. 4. 

The vectors V in Fig. 3 turn 2ir in the direction along a circle around the critical point 
(/q,(?q)- Since the right hand sides of (|6.1|) are regular in t 6 the components 

of V are continuous functions with respect to (f ,go). Then there definitely exists a 
critical point on which V = inside the circle. Therefore this figure strongly suggests 
the existence of the solution with both QK-bolt and T-bolt. 

As commented in section 5, the numerical metric with k = 1 extends over the QK 
manifold M with an orbifold singularity at t = 0, and extends smoothly over the twistor 
space Z at the other endpoint t = ^- If we allow the cone- like singularity at T-bolt, k 
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0.0 1.0 2.0 3.0 4.0 t 



Figure 4: Evolution of a, b and / in t for the critical point (/o>5o) m n = 1 case - 

can take an arbitrary value. The case is commented shortly in Appendix B. 

In table 1 we list numerical values /q and g$ for compact solutions of various dimen- 
sions 4n + 4 of the manifold. As n is increased, f(t) stays almost constant around \Jn + 2 
over the range of t, and the value of &o at ^2 and the interval of the numerical solutions 
t 2 — ti converges. Though the scale of the base space / diverges in the order of ^Jn it 
seems that the 4-dimensional fiber metric converges. It should be noted that during the 
evolutions ///is very small for large n but the terms 4n/// in ()6.1|) contribute in the 
same order oi a/a and b/b. 

7 Summary and Discussion 

In the present work, we have made an investigation of higher dimensional compact Ein- 
stein manifolds. We can view such manifolds as the union of principal >SO(3)-bundles 
over Quaternionic Kahler manifolds. Globally, our compact manifold M is considered as 
a fiber bundle associated with a principal G-bundle P, i.e., M = P x G F (G = SO(3) or 
SU(2)) . The fiber F is a 4-dimensional manifold (orbifold) with the Bianchi IX metric 
on which G acts with cohomogeneity one, and the base space M is the Quaternionic 
Kahler manifold. 

Total space M can be regarded as an evolution of P over a finite 'time' segment 
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n 


Vn + 2 


•J 


ft 




9o 




fo 




bo 


/, 


3 _ fo 


t 


2-*l 


1 


1.73205 


1. 


.75334 





.00253 


1. 


,36425 





.85289 


0. 


38910 


4, 


.65747 


2 


2.00000 


2. 


.05472 





.00441 


1. 


,74912 





.97181 


0. 


.30559 


4 


.49746 


3 


2.23607 


2. 


,30830 





.00633 


2. 


,03541 


1 


.04715 


0. 


,27290 


4 


.37347 


4 


2.44949 


2. 


,52917 





.00811 


2. 


,27774 


1 


.10000 


0. 


,25143 


4 


.28367 


5 


2.64575 


2. 


,72822 





.00968 


2. 


,49323 


1 


.13927 


0. 


,23499 


4 


.21713 


6 


2.82843 


2. 


,91144 





.01107 


2. 


,68983 


1 


.16969 


0. 


,22161 


4 


.16625 


14 


4.00000 


4. 


,07213 





.01767 


3. 


,91020 


1 


.28394 


0. 


,16193 


3 


.98370 


23 


5.00000 


5. 


,06126 





.02091 


4. 


,93003 


1 


.32869 


0. 


,13122 


3 


.91686 


47 


7.00000 


7. 


,04587 





.02420 


6. 


,95105 


1 


.36954 


0. 


,09482 


3 


.85836 


98 


10.00000 


10. 


,03288 





.02616 


9. 


96609 


1 


.39204 


0. 


,06679 


3 


.82719 


223 


15.00000 


15. 


,02219 





.02726 


14 


.97751 


1 


.40431 


0. 


,04468 


3 


.81052 


898 


30.00000 


30. 


,01117 





.02795 


29 


.98879 


1 


.41173 


0. 


,02238 


3 


.80052 


2498 


50.00000 


50. 


,00671 





.02810 


49. 


.99328 


1 


.41333 


0. 


,01344 


3 


.79839 


9998 


100.00000 


100 


.00336 





.02816 


99 


.99664 


1 


.41399 


0. 


,00672 


3. 


.79749 


999998 


1000.00000 


1000 


.00034 





.02818 


999 


.99966 


1 


.41421 


0. 


,00067 


3, 


.79719 



Table 1: Numerical parameters for (4n + 4)-dimensional solutions. Parameters 
/g,<?g for compact solutions, and parameters fo,bo, the difference /q — fo, and the 
interval t 2 — t\ are listed. The cosmological constant is normalized as unity. 



[ti,/^]- If we require the compactness, singularities at the boundaries t\ and t 2 should 
be resolved by nut, Quaternionic Kahler(QK)-bolt or Twistor(T)-bolt. We found new 
solutions numerically which connect QK-bolt and T-bolt in the present paper, where 
T-bolt is characterized by a number k. To the extent of the three unknown variables 
a,b = c, / for the metric form ([2.5)1 . we could complete Fig.2 by using the new solutions 
(broken line) and already known solutions (solid lines). 

Since M has the bundle structure with the fiber F on the base space M then the 
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Euler number is factorized as 



X (M) = x(M) X (F) . (7.1) 
By using the Gauss-Bonnet theorem, x{F) is calculated as 

x( F ) = oTT - 2 / £abcd V ab A n cd = N G (l/2 + 2k) , (7.2) 

OZ7T J p 

where N G = 1 if G = 50(3) and N G = 2 if G = SU{2). The factors 1/2 and 2k in (JZ2D 
represent the contribution from QK and T-bolts, respectively. 

Let us consider the large n limit. Though the scale of the base space / diverges in 
the order of y/n the coefficients ()6.7|) of the local expansion near QK-bolt become 

a 3 = -2g* - A/6 , b 3 = g* - A/6 (7.3) 

and the ones (J6.11j) near T-bolt become 

which might imply that the fiber metric converges to the 4-dimensional biaxial Bianchi 
IX metric 

3oo = dt 2 + a{t) 2 al + b{t) 2 (a 2 2 + of) (7.5) 

in the limit n —>■ oo. 

When the total space connects QK-bolt and T-bolt, the fiber space with the met- 
ric (J7.5)) connects nut and bolt in 4-dimensions. It is known that the 4-dimensional 
biaxial Bianchi IX Einstein metrics connecting nut and bolt singularities is either the 
self-dual Taub-NUT-de Sitter metric or the self-dual Eguchi-Hanson-de Sitter metric, 
and their common boundary represents the Fubini-Study metric on CP(2)[19 . A ques- 
tion is whether the fiber metric coincides with the Bianchi IX Einstein metric in the 
limit n — > oo. Indeed the expansion ()7.3|) locally reproduces the behavior of the self-dual 
Eguchi-Hanson-de Sitter metric near the nut singularity, but the global metric is different 



22 



since the value k = 1 at the other endpoint (bolt singularity) for the fiber metric is not 
allowed for the 4- dimensional solution. 

This situation can be compared with the metric on HP(n + 1) (j HP(n + 1) solved 
by Page and PopejH] in the ansatz (Rl) ; in the large n limit, / stays constant and 
the fiber metric gives the standard S 4 metric i.e., the full metric approaches the direct 
product Einstein metric. In contrast, the numerical metric presented here is not the 
direct product metric even in the large n limit, where the existence of the gauge field A 1 
is important. 

Finally we discuss the relation to the Spin(7) holonomy metrics. In the beautiful 
papers [20J [2 lj Hitchin constructed a family of 4-dimensional self-dual Einstein metrics 
with positive Ricci curvature in a triaxial Bianchi IX form parameterized by an integer 
N. These solutions connect the two bolt singularities characterized by k = 2 and by 
1/N. The metrics are explicitly given by the solution to the Painleve VI equation[23 
and approach the Atiyah-Hitchin hyperkahler metric[22] with Sp(l) holonomy in the 
limit N — > oo. 

In 8-dimensional case, the asymptotically locally conical (ALC) Spin(7) metrics would 
be counterparts of the Atiyah-Hitchin metric [T3* | l2*4* | 12*3] . The local behavior (|4.9|) at T- 
bolt with k = 2 reproduces the ALC Spin(7) metric when one of the free parameters 
is adjusted suitably. Furthermore the expansion (J4.1U)) of the fiber metric at the other 
T-bolt with k = 1/N has the same form as the Hitchin metrics. It is tempting to expect 
that there is a series of Einstein metrics with positive Ricci curvature in 8-dimensions 
which approach the ALC Spin(7) metric in a suitable large N limit. Therefore, it is 
interesting to consider solutions in the general metric form (J2.5)) which connect two T- 
bolt singularities, k = 2 and 1/N. We leave this issue for further research. 

Acknowledgements 

We would like to thank H.Kanno for helpful discussions. This work is supported in 
part by the Grant-in-Aid for Scientific Research No. 14540275 and No. 14570073. 



23 



Appendix A 



In this appendix we calculate the Ricci tensor of the off-diagonal metric (|2.4jl . When we 
use the basis <f> 1 {i — 1 ~ 3) for the fiber metric and the orthonormal basis e° , e a (a = 
1 ~ m) defined by ()2.6|) , the non-zero components are 

4 = -E^- (e^+^t) 

k \ I / k 

+ 4° (3) + 474 £ £^0^0^ > 

k,£,rn,n a, 7 

£tf = " f 7 + (m - 1)^ + E ^7 I (A ' 1} 



/ v 7 2 V l f 
1 



+ ~ 771 E E b ij0kiFa 7 0ijF^ , 

ijki 7 



k 

3 ' 



-Ria - 77^ EE^' ^^^/ 3 • 



Here if/ = (l/2)&ifc& fcj and Rf„° denotes the Ricci tensor on 50(3) ; 



K 3 



SO(3) 



E r ^ r ifc ' r y ~ 2 ( £ijk + E e ^ mbm ^ M + E e 3 f ^rnib M j • (A. 2) 



When we impose the assumption of 4n- dimensional QK manifolds (see section 3) , then 
R ia = and the explicit dependence of the group element 0jj disappears from the 
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equations, 

E E hkbje0mkF™0 ne F^ = An E b ik b kj 

k,£,m,n a, 7 k 

E E M^iA** = E 6 ^ > 

7 i 

which yields the 50(3) invariance of the Einstein equation. 



(A.3) 



Appendix B 

Putting y = b — c, and then using the Einstein equation (|3.1()j) we obtain 

y = py + qy, (B.i) 

where the coefficients p and q are given by 



a a f 

P= 4n- , 

a / 

+ -^(& 2 + &c + c 2 )-A. 



(B.2) 



Taking account of the expansion (J4.8|) . we approximate the equation (|B.1|) in the limit 
t -> 0, 



?/ + 72/-^?/ = ' ( R3 ) 



1 4N 2 

which has the regular solution y = t 2N for a positive integer N. Thus we conclude that 
the expansion takes the form of (J4.1(J|) . 



Appendix C 

Let us consider k ^ 1 case, where cone- like singularity appears at T-bolt. The angle 
around T-bolt is 2irk then the singularity is angular deficit for k < 1 and angular excess 
for k > 1. In this case, solutions make a family parameterized by k. A curve on the 
(fo, 5'o)-plane depicted in Fig. 5 shows a family of solutions. 
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0. 0026 
0. 0022 
0. 0018 
0. 0014 
0. 0010 
0. 0006 

1.5 1.6 1.7 1.8 1.9 2.0 

h 

Figure 5: The curve on (f , g )-plane which denotes a family of solutions with a 
cone-like singularity in n = 1 case. The open circle on the curve, critical point 
(/q, <?q), is the solution with k = 1, which is regular at T-bolt. The left branch with 
respect to the critical point corresponds to the solutions with angular deficit and 
the right branch does angular excess. 
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